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Abstract 



We evaluate matching conditions for the Wilson coefficients of operators mediating the 



Oh! 
> 
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^ 6 — > S7 and b ^ s gluon transitions in a large class of extensions of the Standard Model. The 
calculation is performed at the leading order in flavour- changing couplings and includes two- 
loop QCD corrections. These corrections can be numerically important when approximate 
cancellations occur among the new physics contributions and/or the SM one. 
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1. Introduction 

The decay B—*Xs'~f is known to be a sensitive test of various new physics scenarios. 
Its branching ratio can significantly deviate from the Standard Model prediction in the 
Supersymmetric Standard Model (SSM) |]T|, multi-Higgs-doublet models 0, ^, left- 
right-symmetric models [0, ^, and other theories f^. Since the experimental results of 
CLEO 

BR[B ^ X,7] = (3.15 ± 0.35stat ± 0.32,^,^ ± 0.26rnodei) x 10"^ (1) 
and ALEPH 



BR[B ^ X,7] = (3.11 ± OMstat ± 0.72,^,^ x 10^" (2) 

are consistent with the SM expectation (see below), parameter spaces of these theories 
get severely constrained. However, given the present sizable experimental errors, it is still 
conceivable that new physics effects are large, but either tend to cancel in the decay rate 
or tend to reverse the sign of the amplitude, when compared to the SM. In such cases, 
next-to-leading QCD corrections to the new physics contributions can be numerically very 
important H, gj. 



After the next-to-leading QCD corrections have been found [0-|]T4| and confirmed [|T5 



TTj in the Standard Model case, all the SM predictions given in the literature fall in the 
range 

BR[B Xs-f] = (3.3 ± 0.3) x 10"^ (3) 



All the authors agree that the present theoretical uncertainty is around 10% (see e.g. [14 
jl7|-|^). Even if all the authors had included the same set of corrections (which was not 
the case), the central values of their predictions would be allowed to differ by up to around 
7%, which is the expected size of the uncalculated next-to-next-to- leading QCD corrections. 
Apart from the QCD corrections, two other important contributions have been calculated 



recently, namely the A/m^ non-perturbative corrections ||2T[ and the leading electroweak cor- 



rections [|1^, |2^. In both cases, the effects on the branching ratio are below the ~10% overall 
theoretical uncertainty, but the shifts in the central value are relevant. 



Theoretical analyses of 6 — > 57 and b ^ s gluon are usually performed in three steps. 
First, the full Standard Model (or some of its extensions) is perturbatively matched on 
an effective theory containing only light degrees of freedom, i.e. particles much lighter 
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than the electroweak gauge bosons. Flavour-changing interactions in the low energy theory 
are mediated by effective operators of dimension higher than four. In the second step, 
Wilson coefficients of these operators are evolved with the help of the renormalization group 
equations from the electroweak scale /io ~ to the low-energy scale /if, ~ rrif,. Finally, one 
evaluates matrix elements of the effective operators between the physical states of interest. 

This procedure allows a resummation of large QCD logarithms such as [a^ ln(M^/m^)]"' 
from all orders of the perturbation series. Moreover, its important advantage lies in that 
both the RGE evolution and the calculation of matrix elements are practically the same in 
the Standard Model and in many of its extensions. Thus, a variety of new physics effects in 
B-^Xg'y can be analysed with next-to-leading accuracy without having to repeat the involved 
calculations of three- loop anomalous dimensions [0, two- loop matrix elements [|13| or 



Bremsstrahlung corrections [[ill, |T5|. It is sufficient to perform the matching calculation 



for 6 — > S7 and b ^ s gluon, including the potentially relevant next-to-leading two-loop 
contributions. 

Performing matching calculations in theories containing exotic particles requires calculat- 
ing basically the same (or at least very similar) sets of Feynman diagrams as in the SM case. 
The variety of final results is to a large extent due only to differences in electromagnetic 
charges and QCD-transformation properties of the heavy particles that are being decoupled. 
Consequently, it seems reasonable to first calculate the matching conditions for a relatively 
generic extension of the Standard Model. Results for various specific theories can then be 
obtained by only substituting particular values of couplings, charges and colour factors to the 
"generic" formulae. Presenting such "generic" formulae for the one- and two-loop matching 
conditions is the main purpose of the present paper. 

The class of models we shall be interested in are theories in which the leading contributions 
to 6 — > S7 and b ^ s gluon originate from (heavy fermion)-(heavy boson) loops, similarly 
to the SM case. Each such loop gives an additive contribution to the considered Wilson 
coefficient. Thus, it is sufficient to perform a calculation with only a single heavy fermion 
and a single heavy boson (scalar or vector). Having calculated the one-loop diagrams and 
two-loop gluonic corrections to them, we shall be able to reproduce the known next-to- 
leading matching results in the Standard Model,0 the Two-Higgs-Doublet Model (2IIDM), 
and gluonic parts of the NLO corrections to chargino contributions in the SSM. Examples of 
new results that can be obtained from our formulae are the NLO matching contributions in 



except for the (light quark)-(VF-boson) contribution, which we calculate separately. 



the left-right-symmetric models, as well as gluonic parts of the NLO corrections to neutralino 
and gluino contributions in the SSM. 

Completing the full NLO calculation in the SSM requires, in addition, evaluating two-loop 
diagrams containing no gluons but only heavy particles in internal lines. Some of those results 
are collected in appendix C. The remaining contributions will be treated approximately in 
the heavy gluino case. 

Our paper is organized as follows. In the next section, we evaluate leading and next- 
to-leading matching conditions originating from (heavy fermion)-(heavy scalar) loops. In 
section 3, a similar calculation is performed for (heavy fermion)-(heavy vector boson) loops. 
Light-quark contributions are discussed in section 4. In section 5, we specify what substitu- 
tions have to be made in the general results, for particular extensions of the SM. Section 6 is 
devoted to presenting a numerical example of the NLO matching contribution effects. There, 
we consider the SSM with decoupled gluino, which is relevant, for instance, in the context 
of models with gauge-mediated SUSY breaking. Appendix A summarizes several useful for- 
mulae which we have used for evaluating colour factors in the generic case. Appendix B 
contains a list of functions that enter our "generic" results. Finally, appendix C contains 
results for matching contributions originating from quartic squark vertices in the SSM. 

2. Heavy fermion — heavy scalar loops 

In the present section, we shall evaluate one- and two-loop matching contributions for 
b —>■ S'j and 6 — >• s gluon in a theory described by the following lagrangian: 



^ = CQCDxQED{u,d,s,c,b) + {D^(j))\Dij,(f)) - ml(f)'^(j) + tpli^ - m^)tl; 

+ [{SlPl + SnPR)sk + (BlPl + BuPr^] + h.c.} + (4) 



Here, CQCDxQED{u,d, s,c,b) denotes kinetic terms for the light quarks, photons and gluons 
as well as their gauge interactions. The remainder of the first line contains kinetic, gauge 
interaction and mass terms for a heavy complex scalar and a heavy Dirac fermion. The 
covariant derivatives of these fields are the following: 



The couplings Sl^r and in eq. (^) parametrize Yukawa interactions of the heavy 
particles with the s- and 6-quarks {Pl,r = (1 =F75)/2). These interaction terms must be 



'd, + ig^Gl%-^+ieQ^A,] 0. 



(5) 



3 



QED- and QCD-singlets. The heavy particles can reside in any representation of SU{3)coiour 
for which such singlets exist. The following gauge-invariance constraints must be satisfied 
by the electric charges and colour generators: 

+ Q<f, = 

'^{■^)jnCink + T^,j,)in^njk = CijnT^k^ (6) 

where Cijk are the Clebsch-Gordan coefficients contracting colour indices in eq. (^), and T*^ 
on the r.h.s. above is the generator of the fundamental representation of SU{3). All the 
generators satisfy the standard commutation relations 



rpa rpb 



-la rph 



ifabcTl^^)^ T('^),T(^) — ifabcTl"^), T"',T —ifabcT 



la rpb 



(7) 

The Casimir eigenvalues for the heavy particles will be denoted by and k^, respectively: 

^ '^'A-'-' ^(0)^(^) = (8) 

We assume that the Yukawa interactions Sl,r and Bl r are weak, i.e. that it is mandatory 
to calculate at the leading order in these interactions and aem, but at the next-to-leading or- 
der in as- Furthermore, we assume that all the remaining interactions (denoted by Cirreievant) 
do not influence the b ^ s'j and b ^ s gluon amplitudes at the considered order in pertur- 
bation theory. 

In each particular model, one needs to carefully verify whether the latter assumption does 
not lead to neglecting important contributions. Even in the Standard Model case, the top- 
quark Yukawa coupling Yt and the quartic Higgs coupling A can be similar in magnitude to 
the QCD gauge coupling. Potentially, the 0{Yt, A) corrections to the considered amplitudes 
could be almost as important numerically as the 0{as) matching corrections. However, it 
has been explicitly checked in ref. that this is not the case for b 57, i.e. that 0{Yt, A) 
effects in the SM are only around 1% in the branching ratio, which is partly due to accidental 
cancellations. 

In the present and in the following three sections of this article, we shall neglect Cirreievant, 
i.e. we shall restrict ourselves only to the leading one-loop diagrams, and to the next-to- 
leading two-loop diagrams with gluons^ Such diagrams involving (heavy fermion)-(heavy 

scalar) loops are presented in figs. 1 and 2, respectively. Small circles denote places where the 

In section 6, where the SSM case wih be considered numerically, we shall include also flavour-conserving 
gluino couplings and the quartic squark couplings, which are proportional to ag. 



external gluon can couple in the 6 — > s gluon case. The external photon in the b ^ case 
can couple in the same places, except the ones denoted by "5" on internal gluon propagators. 
Thus, we have 45 two-loop diagrams in the b —>■ s gluon case, and 37 two-loop diagrams in 
the b ^ case. 

^ © 

Figure 1: One- loop diagrams 




Figure 2: Two-loop diagrams 

We evaluate all the diagrams off shell, in the 't Hooft-Feynman version of the background 
field gauge for QCD. For all the quantities we need to renormalize, we use the MS scheme 
with the renormalization scale /Iq that is assumed to be of the same order as the heavy 
masses. Before performing the momentum integration, we expand the integrands up to 
second order in (external momenta)/ (heavy masses). All the spurious IR divergences arising 
in this procedure are regularized dimensionally. 
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Next, we require equality of our result to the similar off-shell IPI Green function in the 
effective theory described by the following effective lagrangian 



where 



4:Gf 

^QCDxQED 

{u,d,s,c,b) + —^V^lVtb 



1=4,7,8 



EOM- vanishing 
operators 



^m,{m^^^PRh)F^ 



Pa 
Pj 

P8= j^m,{sa>^'^T-Pnb)G%, 



(9) 



(10) 



and the primed operators P/ are obtained from the ones above by changing Pl,r to Pr^l- 
The structure of the EOM- vanishing operators (i.e. operators that vanish by the QCDxQED 
equations of motion) as well as other elements of our matching procedure can be found in 
ref. ||2^ where the Standard Model case is described in detail.^ 
The obtained Wilson coefficients 



CM = C^'\^^o) + (/xo) + ... {^ = 4,7,8) 



(as well as their primed counterparts) can be written as linear combinations of various 



functions of a; = (m^/ 



\2. 





= 0, 




= N{Ri 
+R2 




= N{R^ 
+R2 




= NSIBl 




= 3 In 4- 
+N{R, 
+R2 




= 3 In 4 

'Til 

+N{Ri 
+R2 



} 



g2{x)Qi, + (xg2{x) + \ 

- (giix) + \g^ (i)) + \g^{x) - ^gi (i) 
|(k^ - K^) ((1 - x)g2ix) - I) + ^g2ix) + i] } , 



16^(0) 
27 '-'S 



lhi{x)K,^ + h2{x)K^ + h3{x))Q^ + {h4{x)K^+ /i5(x)k<^+ hQ{x)) 
[(/i7(x)/t^ + hs{x)K,^ + hg{x)) + {hio{x)K^ + hn{x)K^ + hi2{x)) Q^]} 



{2k^ - K^)x 



(0) 



dx 

/ii3(x)k^ + hi4{x)K,l + hi^{x)K^K^ + hi(i{x)K^ + hn{x)K^ + hi^{x) 

/ll9(x)/t^ + h2o{x)K^^ + h2l{x)K^K^ + h22{x)K^ + /l23 (a^) + h2A{x) | 



S,l 



{0)lr 



Four-fcrmion operators generated by the photon exchange (such as {sj^PLb){e'~f^e)) are ignored here 
but included in ref. p4[. They are irrelevant to 6 — > 37 and b ^ s gluon, but they matter, for instance, in 



6 



with 



Ri = SIBl + —S^Br and = SIBr"^. 



The primed Wilson coefficients can be obtained from the ones above by simply interchanging 
the chirality subscripts L ^ R. Explicit expressions for the functions fi{x), f2{x), gi{x), 
g2{x) and hi{x)-h24{x) are given in appendix B. 

The symbols cf^^'^ used in the latter two equations denote those parts of the leading- 
order coefficients Cf'^ that are proportional to S^Br. The global normalization constant 
equals to 

N = ^ (19) 

where ^ parametrizes the contraction of the Clebsch-Gordan coefficients: 

CijkCijn = ^Skn- (20) 

The subscripts "S*, 5" at the derivatives of C^^^ and C^^ mean that the products S^Bl, 
SrBr and S^Br^ need to be treated as independent of x when these derivatives are taken. 

It is remarkable that colour factors of the two-loop diagrams could have been reduced to 
quite a simple form, even though we have not specified the representations of SU{3)coiour 
in which the heavy particles reside. Several identities that are useful in performing this 
reduction are summarized in appendix A. 



3. Heavy fermion — heavy vector boson loops 

Here, we are going to consider only QCD-singlet vector bosons V"^, and heavy fermions 
ip in the fundamental representation of SU{3). The lagrangian is now assumed to have the 
following form: 

^ = CQao.QEoiu, d, s, c, b) - (D^V^nD^V) + m^yV;V'' + ieuoV;V,F^^ + ^^(z^ - m^)^ 
+ [v;i^-i^[{aLPL + (yRPR)s + {l3LPL + l3RPR)h] + h.c.} + C,^^^,^^^„,, (21) 

where F^*^ is the electromagnetic field-strength tensor, and 

D^V^ = {d^ + ieQyA^)V, with Qv = -Q^-\. (22) 

The constants (Jl,r and I3l,r are some arbitrary weak coupling constants. If we substitute in 
eq. {^V ^W,i,^t,Q^-^-l,uo^ -2, aR.jSR ^ 0, ^ ~^Vts and {Sl ^ -^V^t6, 
we obtain all the Standard Model interactions that are relevant to the (W^-boson)-(top 
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quark) loop contributions to 6 ^ S7 in the 't Hooft-Feynman version of the background 
field gauge. 

The present calculation differs from the Standard Model one mainly by that it allows 
right-handed couplings of the heavy vector boson to fermions. Such couplings occur for the 
VT-boson in the left-right-symmetric models. Their small magnitude can be compensated 
by the large ratio mt/rrib in contributions to B^X^'y p. 

The diagrams we need to consider now can be obtained from the ones presented in figs. 1 
and 2 by removing all the diagrams with gluon-scalar couplings and replacing the scalar by 
the vector boson. In this way, we obtain 16 two-loop diagrams in both the b ^ s gluon and 
b ^ S'j cases. 

The effective lagrangian with which our full theory is being matched remains the same as 
in eq. (|^). The contributions to the Wilson coefficients we obtain in the present case read 

C^f (yuo) = 0, 
= N 



(23) 



C}°Vo) 

C^°Vo) 
Cf^(/xo) 



N 



crlPi H -^*r(3r {ji{x)Q^ + j2{x)Qv + j3{x)uj) 

rUb J 

Wn"^ {j4{x){Q^ + Qv) + MxP] 
rUb 



rrib / rUb 
-12a;2 + 18x -4 , -25x2 + 29a; + 2 



In 



/^o 



3(1 -x)4 



Inx + 



9(1 -x) 



(24) 
(25) 
(26) 



ix 



dci_ 

dx 



3 ^ 9 ^ 



+N 



o-lPl H cr|j/3ij {ki{x)Q^ + k2{x)Qv + hix)uj) 

nib 



+o-1(3r— {ki{x)Q^ + k5{x)Qv + ke{x)uj] 
nib 



(27) 



In 



j4_ 

ml. 



da 



(0)' 



ix 



+N 



dx 



(0) 



nib 



(^If^L + —(JrISr ) kj{x) + alf3R^ks{x 



mb 



{2t 



As before, contributions to the primed coefficients can be obtained by simply interchanging 
the L and R subscripts in the couplings a and (3. The normalization constant N reads 

V2 



N 



Sml^CpVtlVtb' 



(29) 



The functions ji{x)-j5{x) and ki{x)-ks{x) of x = {m^/mvY are given explicitly in appendix 



B. 
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4. Contributions from (11/^-boson)— (light quark) loops 

In order to obtain complete matching results, we need to consider contributions from 
loops containing light quarks. In the Standard Model case, loops with 14^-bosons and charm- 
quarks are relevant. Diagrams with up-quarks are analogous, but numerically less important 
because of the smallness of the CKM matrix element Vub- 

In the following, we shall assume that only the Standard Model contributions are non- 
negligible in the light-quark case. This is a correct assumption in many extensions of the 
Standard Model (e.g. 2HDM or SSM) in which couplings of heavy scalars to light fermions 
are suppressed by small Yukawa couplings. In the left-right-symmetric models, there are 
additional contributions from right-handed coupling of the W^-boson to quarks. However, 
for the light quarks, such contributions to 6 — S7 and b s gluon are not enhanced by the 
large ratio mt/mf,. Consequently, they are negligible when compared to the top-quark one.0 

The Standard Model diagrams we need to consider here are the same as in the previous 
section, with the heavy fermion replaced by the light quark. However, the effective theory 
side is now somewhat more complicated, since we need to include operators that contain 
light quarks. In the charm-quark case, they read 

P2= {s-f^PLc){c-f^PLb). 

Their Wilson coefficients are found by considering the shcc Green function^ 
Cr(/io) = 0, Cf (/io) = r,„ 



-15 - 61n^ 



C«(/io) = 0. 



(31) 



where = KlKt/V^^l^fc- The value of C(^) has been obtained in the MS scheme (applied 
throughout this paper) and using the evanescent operator 

El = {s^^^.^pPlT''c){cYYYPlT%) - I6P1. (32) 

After performing the matching for Pi and P2, contributions to the Wilson coefficients of 
P4, Py and Pg can be found. We obtain 

5^CP ifio) = 0, H/io) = (1 - fln^) (33) 



^ Unless the left- and right-handed CKM matrices have a very different hierarchy of their elements, which 
we assume not to be the case here. 
^ See ref. IE4| for more details. 



Vo) = %r.u 5^C^J'Vo) = (-1 - ^ In ^) r.t. (34) 

5'^Cf\ii,) = ir,„ S'^ci'\f^o) = {-§-, + ^ In 4;) r^t- (35) 

Results for the up-quark loops are obtained from the above ones by replacing by 
Tut = V*^Vub/V^*Vtb. After adding the up- and charm-quark contributions, we can make 
use of the equality r^j + = —1, which follows from unitarity of the CKM matrix. 

The primed coefficients are negligible in the Standard Model case, because they are pro- 
portional to small Yukawa couplings of the would-be Goldstone boson to light quarks. 



5. Substitutions 

Let us first use our generic results to reproduce the known SM matching conditions at the 
next-to-leading order in QCD. The only scalar we need to consider in the SM is the charged 
would-be Goldstone boson. We use the 't Hooft-Feynman version of the background field 
gauge, in which the (14^-boson)-photon-scalar vertices are absent. The only heavy fermion 
is the top quark, for which one finds 



^ (Slij. . ^SUB^) . I (l ^ ^) . I wUh . ^ (36) 

NalPL = I and = = 0. (38) 

The remaining substitutions one needs to make are — >■ m^, m^,my — ^ Myy, 

|, Q^, Qv —1, uj —2, ft^v I ^nd 0. Adding up the results from sections 2, 3 

and 4, one easily finds 



3 r 1 -'-'I V-^; 6(l-x)4 ^^^-^ ' 36(l-x) 

ci%o) = -7^xHI'''\x) with { 

f7l8J('^^ _ X 1 _|_ -x^+5x+2 



2 



m'i(rr\ — -3x^+2x 1^ ^ I -8x-'-5x+7 
y-^) — 6(l-x)4 n^-^ 36(l-x)3 

[8]/ N ^ , . _^2j 



,^{1)/,, \ _ 4x^-]fe3-|-15^ 1 „ I 7x3-35x^+42a:-8 , 2 Mq 

(--4 l/^Oj - 6(l-x)4 12(l-a;)3 ^ 3^^^' 

^((1)/,, ^, -16g-^-122a:3+80a:^-8x r A 1^ i -4x'^+407x'*+1373x3-957a:^+45a; i^ ^ 

'-7 l/^oj - 9(TZ^)4 L«2(^i--J + 81(l-x)5 



1520x4+1296l2:='-12126x2+3409a;-580 



486{l-a;)4 



6x4+46x3-28x2 



In-^^ 

mf 



1 I 106a:4+287x'^+1230x2-1207x+232 

3(l-a;)5 '"^-^ 81{l-a;)4 

^(1)/,, N _ -4x"t+40x3+41x2+x r„- f^ l\ I 32x5-16x'*-2857x3-3981x2-90x i^ ^ 
<-8 l/ioj - 6(TZ^)4 i.22(^i--J + 216(l-x)5 X 

_l_ 611x-4-13346x3-29595x2+1510x-652 , [ -17x3-31x2 In 'r -I- 89x-^-446x3 - 1437x2 -950x+152 1 i 

1296(l-x)4 [ 2{l-x)5 -\- 108(l-x)'* \ rn^' 

These results agree with the main findings of refs. |T2|, |T6|, |1^]. Contributions to the primed 
coefficients can be neglected in the SM, because they are suppressed by m^/mfc. 
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Let us now turn to the Two-Higgs-Doublet Model II. For diagrams with the physical 
charged higgson exchanges, the necessary substitutions in the formulae of section 2 are the 



following: rut, Mh-, |, Q<i 



■1, — i> |, — i> and 



nsibr 



nib 



y 

2 

y 

2 



cot^ /3 1 



fcot^A 



with y 



Ml. 



(43) 
(44) 



The resulting contributions to the Wilson coefficients read 



with 




3(l-y)3 ^^^y + 6{l-j/)2 



2{l-!/)2 



5«-cf^(^o) 
5«-C?)(/io) 



cot 



+ 



1202^;'^+7569^;3 - 5436^;^ +797</ 



486{l-y)4 

,,3 j_l 1 '),,2_ 



+ 



27(l-y)4 

, -32^/3+112^/^-48;; ^^ fl-il + i%!z:12%lf6% , 8^3-52^^+28^ 



In 4} 



+ 



+ 



~12y3-56y^+32^ 
9{l-y)4 



Iny + 



16^3-94^2+42^ 
9(1-^)3 



In 4, 



cot 



0{ 



13y-*-17y3+30y2 



6(l-y)4 

-4451^^+7650^3_x8i53;^2^^;^]^3Q^ 
1296(l-y)"' 

17^2+19^ 1 , 7y3-16y2+81^ ' 



lV-31^1 , -7y^+18;/3_26i^2_38^ 

6(i„„)5 iU^/ 36(l-y)4 
, , 29;;3_44^2^.i43^ 



y3_3i„2 

6(i-y) 

34y3-7y2^]^g5^ 
12(l-y)4 

In 4. 



In 4) 



These results are in agreement with refs. 0, ^. Similarly to the SM case, the primed coeffi- 
cients are suppressed by ms/mh, and can be neglected. 



The third example we would like to discuss here is the SU{2)l x SU{2)fi x f/(l) model. 
The notation introduced in ref. will be followed. We shall restrict ourselves to such 
contributions from the light PF-boson (and the corresponding would-be Goldstone boson) 
which are suppressed by the small VT-boson mixing angle but simultaneously enhanced 
by the large quark- mass ratio In this case, the substitutions one needs to make in the 
results of sections 2 and 3 are the following: 



(49) 
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where x 



and 



92L »/tf,(^) 



(50) 



The "left-left" and "right-right" products of couplings can be set to zero now, because their 
non-SM parts are of order 0(C^, 1^)' negligibly small. Substitutions for masses, charges 
and colour factors here are the same as for top-quark loops in the Standard Model case. 
The obtained contributions to the Wilson coefficients read 

(51) 
(52) 

(53) 



5""Cf)(/xo) 


= 0, 
= A 


^--CfVo) 


= A 


^""C^j'Vo) 


= A 



tb 



tb 



1„ ™ I -5a:2+31x-20 
LLLX -\- i2(l-x)2 



3x'^-2x 

2(l-a;)3 "^"^ ' 12{l-x) 

—3x 1^ ^ I —x^—x—4 

2{i-xy 



Inx + 



4(l-x)2 



-32x3-112x2+48a; 



+ 



+ 



9{l-a;)3 
24x3 _|_3203^.2 _ 220X+20 
9{l-x)^ 

■nx^+89x'^+12x 



86x3+1202;2-30x-32 



Inx 



+ 



9{l-x)i 

12x3+563;^-32a: 1^2;+ 16x3+90x^+66x-64 



6(l-xp" 

87a:'^ -640x^-451x- 148 
24(l-x)3 



3(l-x)4 

(l - + 

+ 



9(l-a:)3 



1^4}' (54) 



34z3-375x^-207x-28 
12(l-a;)4 



Inx 



17x2-19xj^^ _|_ _7x'3-36x2-159a;-28 



ln4|. 



(55) 



^^^-^ ^ 6(l-x-)3 

The primed coefficients are obtained from the above ones by changing A^^ to (A*'^)*. Expres- 
sions for the leading-order coefficients are in agreement with ref. 0. The next-to-leading 
results are new. 

Physical charged scalars present in the left-right-symmetric models can give contributions 
enhanced by too 0. Such contributions can be calculated with the help of our generic 
formulae as well. However, their explicit form and potential numerical importance depend 
on details of the Higgs sector that is not unique in these models. 



Finally, let us turn to the Supersymmetric Standard Model. In this case, we shall exactly 
follow the notation of ref. In particular, the relevant couplings of charginos, neutrali- 

nos and gluinos to matter will be parametrized by X^-^'^, Z^'^'^ and g^T^^''^, respectively. 
However, we shall formally treat all these matrices as independent from other parameters 
of the model. Only diagrams containing gluons have been included in their one-loop QCD 



renormalization, and we have used the MS scheme in dimensional regularization (not in 
dimensional reduction). In effect, immediate substitution of tree-level supersymmetric ex- 
pressions for these matrices (eq. (2.12) of ref. ||2^) is allowed only for the leading-order terms 
The CKM matrix will be denoted by K in all the SSM formulae. 
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in the results presented below. Section 6 contains a discussion of applicability and treatment 
of the 0{as) terms. 

In the results presented below, flavour-conserving gluino interactions as well as strong 
quartic squark interactions have been formally treated as the weak ones. Additional con- 
tributions to the considered Wilson coefficients from two-loop diagrams with strong quartic 
squark interactions are collected in appendix C. A complete inclusion of strong gluino cou- 
pling effects at NLO is beyond the scope of the present paper. Consequently, our results 
can be used for making numerical predictions in the SSM only in the heavy- gluino limit (see 
section 6). 

At one loop in the SSM, & — > 57 and 6 — > s gluon off-shell Green functions receive 
contributions from chargino-squark, neutralino-squark and gluino-squark loops, in addition 

to the SM and charged higgson contributions that have been already discussed. For the 
loops containing chargino xj (-^ = 1)2) and the up-squark -u^ (A = 1, .... 6), as well as for 
the gluonic corrections to them, we need to make the following substitutions in section 2: 
m^^m-j, m^^ma^, ^ -1, Q^^2/3, ^ 0, and 

After summing up all the chargino and up-squark species, we obtain the following contribu- 
tions to the Wilson coefficients 

92J^tsJ^th A=ll=l 



X 



A3 



A3 



The constants A'^l = | and A'^l = | drop out by unitarity of squark mixing matrices when 
the chargino-matter couplings X^'^'^ are no longer treated as arbitrary but expressed in 
terms of other SSM parameters. The functions Hi2^ have already been encountered in the 
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SM and 2HDM cases. The functions if| , if| 2 ifj 2 have the following explicit form: 

ffWrr^ - -1 In T -^ -2x2+7x-ll 

-ni I^J-; — 3(^_^)4inx ig(i_j.)3 5 

C7[7]'/t,^ — 24x3+52x^-32^ t A ( ^ _ }_\ i "189a:3-783x^+425x+43 i^ ^ , - 1030x3 - 1899x^+1332^+85 

-"l y-^l ~ 9(l-x)4 -^'2 x^ 81(l-x)5 "^-^ 243(l-x)4 ' 

U-I^l'Va^^ — 6x3-62x^+32x 1^ ^ I 28x3-129x^-12x+41 

J^i \X) — 9(i_2,)5 niJ- 27{l-x)4 ' 

[7]'/^^ — 112x2-48xr- i lA i 12x3-176x^+64x+16 i^ ^ , -170x^+66x+20 



IJL'I /^N _ 112x'-4Sx T A M 1 ^ i 12x-^-lYbx"+b4x+lb 1^ ^ 

^2 KX) - -g^Y3^i.Z2 (^i - -j H 9(iz^)4 ma; -r g^j^^p , 

It['']"/'t,^ — 12x3-88x^+40x 1^ ^ I -14x^-54x+32 

-"2 l-^J — 9(l-x)4 "^-^ 9{l-x)3 ' 

fjISlVa-'l — -9x3-46x^-49x r„- A _ i ^ i 81x3+594x^+1270x+71 ^ , 923x3+3042x^+6921x+1210 

-"l \-^) ~ 12(l-x)4 ^''2[i- -r 108(l-x)5 "^-^ 648{l-x)4 ' 

TTi^]"{rr\ — 5x^+19x i„ ™ I 7x3-30x^+141x+26 

-"1 y-^) ~ 3(l-x)5 "^-^ 18(l-x)4 ' 



3(l-x)3 -^-^ I 6(l-x)-* ' 12(l-x)3 ' 

tt[^]"(^\ _ 10x^+26x1^™ I -x^+30x+7 
-"2 l^J - 3(l-x)4 ^^^-^ 3{l-x)3 • 

Contributions to the primed coefficients can be obtained by interchanging X^'^ and in 
eqs. (^) and (^). Their suppression by rris/mb becomes transparent only after expressing 
X'^^ in terms of other SSM parameters (cf. eq. (^) below). 

Very similar substitutions need to be made for loops containing the neutralino {I = 
1,2,3,4) and the down squark {A = 1, ...,6): — > m^o, rriii, — > m^^, —^0, — > 
— 1/3, —^0, 4/3 and 

V ^ ^ m, ""J glKlKt, ^ ' I A2 \ ^ J as J ' 



NSlBn^ = , ' ^ (zF^y (ZF-) ^. (6^ 



One could ask whether the results obtained in section 2 for heavy Dirac fermions can be 
applied to the case of neutralinos, which are Majorana particles. The answer to this question 
is positive: for the particular set of Feynman diagrams we have considered, there is technically 
no difference between Dirac and Majorana fermions on internal lines. No "clashing-arrow" 
propagators have to be included, and no extra combinatoric factors occur. 

Application of the above substitutions gives us the following neutralino contributions to 
the Wilson coefficients: 



14 



5* C^7,8(/io) 



EE 



^ /lf2 



glKlKtb m|o 



X 



A3 



A3 



r[7,8]' 



r[7,8]". 



[7,8]', 



where =i/r(x), i/f (x) = iff + iY; 

= (^), h'^Xx) = -\ {Hf\x) + 1) and 



j4_ 

2 

OA , 



, (70) 



[8]", 



X] 



r[7]' 



9(1 






-28a; i 


9(1- 




16x2 


+48x 


9(1 




8a;2- 


-44a; i 



LZ2 1 



I -108a;2-358x-38 !„ ^ i 23a;^-765x2-693x-77 
81(l-a')^ 243(l-x)4 ' 

r'^-|-42x2-84x-22 
27(l-x)'l ' 



I 10x2-30x-16 



-26a;2-54a;-4 
9(l-x)3 ' 



(71) 
(72) 
(73) 
(74) 



^^4 ^o.; — 9(i_2.)4 ii^-^ 9(l-a;)3 • 

Contributions to the primed coefficients can be obtained by interchanging Z^^ and Z^^ in 



eqs. (69) and m 



Finally, we turn to contributions from gluino-squark loops, i.e. from one-loop diagrams 
with gluinos and two-loop diagrams with both gluinos and gluons. The heavy fermion and 
scalar considered in section 2 are now in the adjoint and fundamental representations of 
SU{3)coiour, respectively. Comments concerning Majorana fermions we have made in the 
context of neutralinos apply here as well. The necessary substitutions in the present case 



are the following: 



m. 



m 



dA ' 



Qt(; ~^ 0, —1/3; S 



N SIBl + —S^Br 
rrib 



NSlBn"^ 
rrih 



89l 



3, — * 4/3 and 

,2 \ 



A3 



+ O 



m2 



89l 



^ 



Ja2\ JASrrib 



The resulting contributions to the Wilson coefficients read 
89l 



2>giKlKt, m 



a A=i 

'9 A=l 



with Va 



m- 



m 



2 ' 



A2 V /A3 



[7,' 



(va) In 




(75) 



(76) 



15 



A3 



47r 



H, 



[7,8]', 



[7,8]" 



Va) In 




where i/f (x) = - i/fp' (x) , (x) = - 1 (iff (x) + |) and 



hF(x) 



i7f(x) 
i/P'(x) 



H, 



[8]'' 



18x'^-27x^+l 1„ ™ _|_ 73x^-134^+37 
24(l-x)4 ^^^-^ 72(l-x)3 ' 



17x^+86: 
18( 



f^L^, (1 - i) 



12{l-x)5 



I -36x^-3153:^+1161x^+751a:+23 i , 
162(l-xO^ ~^ 

18x3+107a:^+43a; 1„ ^ _i_ 
18(l-x)5 ■'^^■^ 



In X 



-799x3+1719x2+10431x-1847 



972(l-x)4 



-5x^+384a:^+609x+20 
108{l-x)4 ' 



19x^+60x-15 
9(l-x)3 



6x^+36x^+48x-5 i„2 
6(l-a;)4 



+ 



(l - ^) + 

27x'^+106x^+52x+l 3; _|_ 14x^+333x-83 



9(l-x)4 ^^^-^ ^ 18(l-x)3 ' 

ISx'^ +80x^+283; ^ , 55x^+69x+2 
9(l-x)4 ^"-^ 9(l-a;)3 ' 



X 1 I 

L6(l-x)4 + 



19x-'+40a:-ll 
16(l-x)4 "^-^ ' 96(l-a;)3 ' 

45x3-1208x^+901x-570 



^LZ2 1 



+ 



96(l-x)'* 

2520a;-* - 10755x'3 - 10638x2 -6427x-44 



-237x3-846x2+282x-95 1^2 



32(l-x)5 
^.3 



In X 



864{l-x)5 



i„ „ _l_ 5359x3-241425x2+143253x-59251 
^^^■^ 5184(l-x)4 



-747x3-640x2+43x !„ ™ 1 



-779x3-7203x2-93x+ll 
288(l-x)'i ' 



9x2— X 1 I 13x— 5 

8(l-x)3 "^-^ 8(l-x)2 ' 



-359x2+339x-204 
24(l-x)3 



(1 - ^) + 



I -207x3-1777x2+23x-151 ^ 1 
48{l-x)4 ^"-^ 



-78x3-333x2+105x-34 1 2 
8(l-x)4 

- 1667x2 +990x- 379 



24(l-x)3 



-126x3-133x2+7x 



Inx + 



-553x2+84x-35 



(77) 



(78) 

79) 
80) 

81) 
82) 
83) 

84) 
85) 
86) 

87) 
88) 



6{l-x)4 ^^^-^ I 12(l-x)3 

Contributions to the primed coefficients can be obtained by interchanging T^'^ and T^^ in 
eqs. (|7|) and (^. 

Our leading-order SSM resuhs agree with refs. |1|, 0, |2^. The next-to-leading results are 
new, except for the chargino ones, which will be compared with those of ref. in the next 
section. 



6. Numerical size of the NLO corrections 

In the present section, we shall give a numerical example of the NLO matching effect in 
the B Xs'~f branching ratio. From among many possible extensions of the SM to which our 
results can be applied, we choose the Supersymmetric Standard Model, because of its current 
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popularity. Our present results allow us to make the NLO prediction for BR[B Xg'y] in 
the SSM only in certain regions of its parameter space. We cannot consider situations where 
the gluino mass is close in size to Mw or nit, because two-loop matching diagrams with no 
gluons (but only gluinos) have not been calculated so far.| In addition, we have to assume 
that tan f3 is not much larger than unity. Otherwise, two-loop diagrams containing no QCD 
interactions at all would become numerically relevant, and the analysis would be much more 
involved. 

For simplicity, we shall consider a scenario in which the gluino is much heavier than all 
the other SSM particles, and we shall neglect all the l/rrig effects. This scenario is somewhat 
different from the one considered in ref. 0, where some of the other superpartners were 
considered heavy as well. 

So long as only the gluino is heavy and decouples, QCD corrections to loops with other 
supersymmetric particles are given by two-loop diagrams with gluons (calculated in the 
previous sections) and by two-loop diagrams with strong quartic squark couphngs (found in 
appendix C). One-loop gluino corrections to the Wilson coefficients of four-quark operators 
vanish in this limit. However, we need to take into account the fact that the theory with 
decoupled gluino is no longer supersymmetric. Consequently, the usual tree-level expressions 
for chargino and neutralino couplings are affected by 0{as) corrections, which do not vanish 
in the rrig oo limit. 

When the gluino mass rrig is much larger than masses of all the other SSM particles, 
one should, in principle, perform the decoupling in two steps. First, one should match the 
complete SSM with the effective theory built out of all the SSM particles except for the 
gluino. This matching should be performed at the renormalization scale fig that should be 
of the same order as rrig. Next, one should perform the ROE evolution of all the couplings in 
the obtained effective theory down to the scale fiQ that should be of the same order as Mw 
or rrit. At this scale, all the remaining SSM particles with masses 0{Mw, rrit) are decoupled, 
and the Wilson coefficients C4, C7 and Cs are found. 

The scales fig and /^o can be set equal to each other so long as In ^ is not much larger 
than unity. This can happen even in situations where the gluino is still heavy enough, for 
instance when rrig ~ 700 GeV, and all the remaining SSM particles have masses below, say, 
350 GeV. The results can then be written in a compact form, because no RGE evolution 

^ It would require including one-loop gluino corrections to the Wilson coefficients of four-quark operators, 
too. 
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occurs between /ig and /iq. We shall take advantage of this opportunity in our example here. 

Once the gluino is decoupled, chargino and neutralino couplings X^^-^, Z^^-^ are related 
to the other parameters of the model as follows (cf. eq. (2.12) of ref. [^): 

Mu 



X 



Ul 



92 



X 



Ur 



-a,v7ir^^ + ay V72r^« 



v2Miy smp_ 



92aYUnT''^K 



V2Mwcos(3' 



92_ 

V2 



N*, + itan^iV;,) + aYN*,T^-—^ 



Mw cos P 



where 



1 + 



92 

V2 



^^(/Uo) 



TT 



■ tan^AT^r^^ + ayNnV^^ 



M, 



D 



M\Y COS (3 



(89) 
(90) 
(91) 
(92) 



In!^ 



12 



and 



ay 



tts(/io) 



TT 



(93) 



In the above equations, all the parameters (including the matrices X^'^'^ and Z^'^'^) are 



assumed to be MS'-renormalized in dimensional regularization, in the effective theory con- 
taining no gluino. Setting a^ and ay to unity, we would obtain relations for parameters of 



the full SSM, MS'-renormalized in dimensional reduction. Explicit values of Ug and ay have 
been obtained by performing simple one-loop matching between these two theories. 

One might wonder whether the appearance of a^ and ay is the only effect of gluino 
decoupling. Potentially, interactions different from those of the SSM could be generated 
at one loop. A more detailed examination of the relevant one-loop diagrams leads to the 
conclusion that other effects indeed occur, but they are never relevant to the NLO QCD 
corrections to BR[B Xs^i]. 

Once ay and a^ from eq. (|93D are included, we are able to verify that our results for 
the NLO corrections to chargino loops (eqs. (|5BD-(|57D, ( |106D and ( |107D ) agree with those 
of ref. [0, so long as the gluino is assumed to be much heavier than all the other SSM 
particles. When performing the comparison, one needs to take into account that our results 



are expressed in terms of MS'-renormalized masses, while the on-shell squark masses were 
used in ref. 0. 

The last worry one could have in the context of the SSM is whether two-loop effects 
involving no couplings proportional to a^, but only to the large Yukawa coupling Yt, could 
be of a numerical importance similar to the NLO QCD corrections. This can potentially 
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happen even in the small tan (3 regime. No definite answer to this question can be given until 
the SM calculation of ref. ||22| is generalized to the SSM case.0 However, the main purpose 
of the present section is only to demonstrate that 0{as) effects can be sizable by themselves, 
independently of the magnitude of two-loop higgson/higgsino contributions. 



Having made all the necessary assumptions, we are now ready to collect the results from 
section 5 and appendix C to test the size of 0{as) corrections to BR[B ^s7] for some 
particular values of the SSM parameters. A relatively simple set of SSM parameters, which 
is allowed by all the experimental constraints (including B — > Xs7), can be chosen as follows: 



tan (3 



3, niyy/ II = —2, 

100 GeV, mgi^ino = 700 GeV, 

140 GeV (the lighter chargino mass). 



{MDab = {Sab - SagSbg) x (350 GeV)^ + SagSbg x (110 GeV)^ 



^AB - Sab + (cos 25" - 1) (5^3(^53 + SagSbg) + sm25° {SagSbs - Sa3Sb6), 



r 



D 



(fiavour- diagonal matrix). 



The first three parameters determine the chargino masses and mixing angles. Down-squark 
masses, neutralino masses and mixing angles are irrelevant here, because we have assumed 
approximate vanishing of fiavour violation in the down squark mass matrix. Consequently, 
all the neutralino contributions to the Wilson coefficients C4.7 § vanish. As far as the relevant 



SM parameters are concerned, we take the same values as in ref. [Tj], i.e. rrit^poie = 175 GeV 



5^A^Cf^(MvK) = -0.197 
6^~C?\Mw) = -0.279 
(5*"Cf^(Mvi/) = 0.272 


6-'^^''Ci°\Mw) = -0.098 
6^'Ci°\Mw) = -0.211 
6^~ci°\Mw) = 0.148 




6^^'^Ci^\Mw) = -2.49 
6^'Ci^\Mw) = 5.24 
6^'Ci^\Mw) = -0.24 
SfC^^^\Mw) = -2.22 


6^^'^Ci^\Mw) = -2.22 
6^'Ci^\Mw) = 2.85 
6^'Ci^\Mw) = 0.18 
6lci'\Mw) = -1.49 


6^^^''dl\Mw) = -0.42 
6^~Ci^\Mw) = 0.02 
6^'dl\Mw) = 0.08 



Table 1. Numerical contributions to the Wilson coefficients in the considered example 



^ The smallness of the net effect in the SM is partly due to accidental cancellations, which may not take 
place in the SSM. 
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and as{Mz) = 0.118. For the Wolfenstein parameters, we choose A = 0.22, A = 0.83, p = 
and rj = 0.3. Such p and t] are roughly in the middle of the allowed region in the SM. They 
are acceptable also in the SSM, for the parameters we have specified above. 

The contributions to the Wilson coefficients obtained at po = M\y are given in table 1. 
The 0{aslnm'^/ pI) terms in chargino and neutralino couplings have been included both in 
the leading and next-to-leading contributions to the Wilson coefficients, as should be done 
in the SSM with decoupled gluino. 

Having found the SSM matching conditions numerically, we evaluate BR[B -^s7] 
according to the formulae of ref. |M and using the same input parameters. However, we 



choose a different photon energy cut-off 6 = 0.90 [|T9l, change aem("^6) to a"^^ c:^ in 
the overall normalization [0, and include the l/rric corrections [^. The final result for 
BR[B X,7] is 

BR[B X,7] = 3.15 X 10-^ (94) 

which is equal to the central value of the CLEO measurement quoted in the introduction. 
On the other hand, if we did not include the two-loop gluonic SUSY contributions (i.e. if we 
set 5^'Ci]l{Mw) and 5^~ci]l{Mw) to zero), we would obtain BR[B X^-f] = 3.69 x 10"^ 
The difference between the two results is close in size to the present experimental uncertainty. 
Such a size of the two-loop SUSY effect is rather generic for light supersymmetric particles, 
for which the leading SUSY contributions to i? ^ Xg'j are of the same order as the SM one. 
Thus, QCD corrections to superpartner loops are expected to be important when making a 
meaningful comparison with experiment in such cases. 

Obviously, satisfying the experimental B —>■ X^'j bound for light superpartners requires 
a certain adjustment of the SSM parameters. In our case, the adjusted parameter was the 
lighter chargino mass. The experimental B —>■ Xg'-f bound would be satisfied within la for 
this mass ranging from 110 to 175 GeV, i.e. no real fine-tuning was necessary. 

For larger tan/?, SUSY contributions to B ^ Xg'y can be much larger than the SM one. 
Then, the NLO SUSY effects are much more important than in the example presented above. 
However, satisfying the experimental B Xg'y constraint requires real fine-tuning then. 

Our results for the two-loop SUSY contributions to B —>■ X^'j could be used in scans over 
the SSM parameter space. The allowed regions in this space would change their position 
for light superpartners. However, our feeling is that making such scans at present would be 
somewhat premature. The NLO SUSY contributions to B ^ Xgj can become qualitatively 
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important only when either hght superpartners are discovered or their existence is almost 
excluded by data combined from many experiments. One of these options will certainly 
be realized when the LHC starts collecting data. At that time, one might appreciate the 
usefulness of the long analytic formulae contained in the present paper. 

7. Summary 

We have calculated matching conditions for operators mediating the 6 ^ 57 and 
b ^ s gluon transitions in a large class of extensions of the Standard Model. Both the 
leading one-loop diagrams and the gluonic corrections to them have been included. Taking 
the Supersymmetric Standard Model as an example, we have checked that QCD corrections 
to new physics contributions can be close in size to the present experimental uncertainty, 
even when the SSM parameters arc not really fine-tuned. A similar situation is expected to 
occur in other theories of new physics containing exotic particles at the electroweak scale. 

The main purpose of our paper was to present complete analytic formulae for the NLO 
Wilson coefficients in a possibly generic extension of the SM. Prom these results, we could 
reproduce the known matching conditions in the Standard Model, the Two-Higgs Doublet 
Model, and for the chargino contributions in the SSM. New results were obtained for the 
left-right-symmetric model, as well as for neutralino and gluino contributions in the SSM. 
Our SSM results form a major contribution to the complete SSM calculation. They become 
complete by themselves in the heavy-gluino limit and for tan (3 of order unity. 
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Appendix A 

In this appendix, we summarize several useful identities involving colour generators and 
the Clebsch-Gordan coefficients Cijk- Thanks to them, all the colour factors in section 2 
could have been expressed only in terms of k^, k^, and ^ (cf. eqs. (|]) and (PPP): 



rjia rjia 



rpa rpa 
rpa rpa 



C* rpa 
ijn-'- (<l))U^ljm 

r-i* rpa r~i 

/ rpa rph \ 

where 



1 [A 
- I - + 



2 V3 
3 

8 
3 

8 



4> I ^ikmi 



3i /4 



3i /4 



f^tp^ab^nm + ~ + '^0 ) f^^bcT^^ + —f]ipdabcTj 



-ic 

■ nm 



(95) 
(96) 
(97) 
(98) 
(99) 
, (100) 
,(101) 



V<t> — 3(k0 - + -K^ - -K^ - -, 

- -^r ^2^3 7 2 



(102) 
(103) 



All these identities can be derived from the basic constraint given in eq. (^). 

For completeness, let us quote the standard identities for the fundamental representation, 
too: 

i 1 1 

(104) 

5 



I, i 1 1 

2/0 



fabcf( 



abd 



36. 



d-abcd-abd g'^cd; 



d. 



aab 



0. 



(105) 



Appendix B 

Here, we present explicit formulae for the functions fi{x), gi{x), hi{x), ji{x) and ki{x) 
introduced in sections 2 and 3. All these expressions are available in Mathematica format 
via anonymous ftp from [ftp: / / feynman.t30.physik.tu-muenchen.de/ pub / preprints/ tum-hep 



321-9 pTunctions.m | . They read: 



fi = STTT^lna; + 



l-x f ' 8(l-x) ' 



h - fin -.^^4 mx + 



6(l-x)4 ^^^-^ ^ 18(l-x)3 ' 
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2(l-x)4 ^^^-^ 12(l-a;)3 ' 

1 In -I- -x+3 

-a:+3 r.- ( T _ l\ I -6a;^+9x^-41a;+6 i„2 , lla:^-32a:3+222a:^+82a:-31 1^ ^ _L 2j;3+lla;^+172a;-17 
3(l-s)2-^*2 1^-1- 12(l-x)5 18(l-a;)B 12(l-a;)4 ' 

4j3+4x2+ll2x+12 7- • /i _ l^ I -6j;3-39x2+71a;+6 i„2 , llj3-150x2-39x+62 !„ „ i -5a;2+76x+481 
9(l_a;)4 ^l'2\^i- ^) -t 9(i-x)5 ^n X -\- 27(l-x)4 m X "h 81(l-x)3 ' 



-Sx r /"l 1^1 51x3 +60x2-1-3; 1^2 ^ I -ll3;'^-163a;3-1533:2-|-353;+4 i^ ^ , 



19x3-163j:2-59a;-|-ll 



12(l-i;) 



4 



7x r.- A _ 1^ I -3a;3-30j;2+a; i 2 , 2x'*'+ll2x^+n7x'^ -Zbx-i !„ ™ i -4a:3+161a2+16a:-5 
6(l-x)2-^'2 l^J- 12(l-a;)S X "h 18(l-x)5 ^ 12(l-a;)4 ' 

-14a;3+112a;2-|-34x r.' A l\ i -39x3+48x2 +23a: 1^2 ^ , -34x3-168x2+78x+8 i^ ^ , -113x2+643x+22 



-112x2+34x T ^ A l\ I -39x3+48x2 +23x 1^2 ^ , -34x3-168x2+78x+8 i^ ^ , 
J3^y4 /.Z2(^i--j + 9(YZ^)5 in X + 27(l-x)4 ^ + 



9(l-x)'4 -^"^ V"^ x) ' 9(l-x)5 ' 27(l-x)'' "^"^ ' 81(l-x)3 

4x-10 TA (^ _ l\ I -6x2-51x-28 i„2 , 6x3-8x2+77x-3 !„ „ i -x2+3x+142 
3(1-^)2 -^^^2 \^-L ^) ^ 6(l-x)4 in X -h 3(i_x)4 -h 6(i_^)3 , 

-2x+8 TA {^ _ }_\ I 6x2-3x+20 i„2 , -6x3+19x2 -85x+9 !„ ™ i 23x2-90x-59 
3(l-x)2 -^*2 -h 6(l-x)4 nl 3(l-x)4 "^-^ ^ 6(l-x)3 ' 

-8x2-8x-176 (^ I 12x2+54x-112 i„2 ^ , -24x2+124x-24 i„ ^ , 34x-104 
9(l-x)3 ^«2(^i--j H 9(1-^)4 in X H 9(i_3;)3 + 9(l-x)2 ' 

-5x-l r.- (^ _ l.\ I -39x2-39x-7 i^2 , 3x3+40x2 +38x-9 !„ ~ i -25x2+147x+22 
3(l-x)2 -^'2 x; 6(l-x)4 3(l-x)4 n^-^ 6(l-x)3 ' 

7.T-1 TA (^ _ l\ I 9x2+9x+5 1^2 I -3x3-29x2-40x+9 !„ ™ i 5x2^138x+7 
3(l-x)2-^*2 \^J- ^) ^ 6{l-x)4 3(l-x)4 niX 6{l-x)» ' 

28x2-176x-44 r,- /i _ l\ i 54x2-66x-34 i 2 , 12x2+100x-36 !„ ™ i 22x-92 
9(l-x)3 ^''2\i- ^) ^ 9(l-x)4 9(l-x)3 "^-^ 9(l-x)2' 

9x-6 r • „ i^ _|_ 45x2+42x-6 i„2 , -3x3-90x2-69x+18 1^ -L 39x2 -183x 
16(l-x)2 -^^2 l^-L ^ 32(l-x)4 32{l-x)4 32(l-x)3 ' 

9x-6 7- • /l _ i\ I -3x2+36x-6 i„2 , 3x3 -45x2 -30x+9 ^ . 3x2-72x+6 

16(l-x)2 -^'2 ^) ^ 32(l-x)4 16(l-x)4 16(l-x)3 ' 

-9x+6 TA (a _ l.\ I -21x2-39x+6 i^2 , -3x3+180x2+129x-36 !„ ™ i -45x2+327x-12 
8(l-x)2-^'2 \^-L ^) ^ 16(l-x)4 32(l-x)4 32(l-x)3 ' 

-172x2+211x-27 



2x2+211x-27 T A f-\ 1^1 Ml^i+639x2+649x-27 i„2 ^ 
96(l-x)3 -t-^2(^i--J H 192(l-x)5 in ^ 



_|_ - 185x^ - l!J9(-i.t.- ' - 1755.t.-- -898X+226 {jix + -50!).t.- '+K5.t.--2(i9!)x-+71 

140x2- 



288(l-x)5 ^^^-^ ^ 192(l-x)4 ' 

-179X+27 r.- /l I -147x3+177x2-569x+27 1^2 ^ 

96(1^^P ^^2(^i--j H 192(l-x)5 ^n ^ 

_l_ 149x4-H072x3+2247x2+790x-226 1^ _|_ 485x3 +479x2+1795x-71 
288(l-x)5 "^-^ 192(l-x)4 ' 

68x3+167x2+314x+75 ta f 1 i\ _i_ -171x3+153x2 +455x+75 ^„2 ^ 
72(l-x)4 -^^2(^1--; H 144(l-x)5 ^ 

_l_ -65x3-381x2+318x+10 1^, _l_ 1513x2 +2782x+1621 
216(l-x)4 "^-^ 1296(l-x)3 ' 

-9 TA {t I -33X-21 1 2 I -3x2+45x+6 i^ „ , -3x+15 

8(l-x)-^^2 l^-L 16(l-x)3 8(l-x)3 "^-^ 2(l-x)2' 

-9 TA A _ i\ _L -3X-15 1 2 , -3x2+45x i , 9x+33 

8(l-x)-^*2 \^-L x; ^ 16(l-x)3 ^ 8(l-x)3 ^ 8(l-x)2 ' 

9 TA A l^ _L 9X-I-9 l„2 ^ , 3x2-45x-3 i„ ^ , 3x-93 



4(T3x)^^2 (l - i) + #^ln^x + '-i^lnx + ^ 



2 ) 



23 



^22 - 24(l-x)^^^2 x) H 12(1-2^)4 ^ 48(l-x)4 ^ 24(1-0;)^ ' 

I, _ -43^+52r- /i l\ I _9i2+37 , 2_ , -63a;=^-2033:^-781a;+39 i^ ^ , -67a;^-2103:-227 

^23 - 24(n^-^^2 - -j + 12(l-a;)4 ^ 48(l-a;)4 ^ 24(1-^)^ ' 

I, _ -7a-^-49rr~112 To 1 _ 1 ^ i 6x^-6x-46 1^2 ^ , 3x-^+70x-39 1^ ^ , -23j:-95 

'^24 — I8(l-a;)3 ^'■2 [i- ^ ) 9(l-x)4 X "h 36(i_a;)3 "1 2; ^ 36(l-a;)2 ' 



• _ -2x^+0: 1„ ™ I -5a:^-5x+4 
-'1 ~ "^-^ 6(l-j;)3 ' 

• _ -x^ In T + -2a;^-5j;+l 
J2 — (i_2:)4iU-i 6(l-a:)3 ' 

J3 = (l-x)3 + 2(L-ji)2 ' 

• _ -4j; Iv, I -2x-2 
J4 — (i_^)3 -h (^_^)2, 

•^5 = (i-x)2 + (T^' 

i, _ 24x3+144x^-80x T ^ A 1^ i 32x3+120x2-56x i„2 

'^i - 3(1^^)3 ^^2(^i--j + 3(iz^)5 in X 

_l_ eOa^ - 1740 j3 _84a:^ +796x- 184 !„ ™ _i_ -1538a:3 -56523:^ +4806j:- 1072 
27(l-a:)5 81(l-i;)4 ' 

7. _ 32a;3+72x^-16x r • l\ , 112x2-16x 1^2 

I 96j;^-1272j3_48a:2+40a;+32 i , - 2024a;3 -954a;^ - 5403:+62 

27(l-a:)B 81(l-a;)4 ' 

7. _ -56x^ r- /i A I -8a;3-64a;^ i^2 , 48^3+80x^-320: i^ „ , 110x^-12a:-2 
'^S - 3(i_^)3^«2 - 3(l-x)4 in ^ 3(l-x)4 H 3(i_2,)3 , 

I, _ 64x^+192a; t A A l\ i 64^!+224x 1^2 ^ , 16x3-368x^-64x+32 i^ ^ , -136x^-256a;+8 
- 3(l-x)a -^^2 - + 3(l-x)'4 ^n ^ 3(r:^P ^{1^ ' 

J, _ 64x2+192xr- /i lA I 64x2+224^ i„2 ^ , 16x3-400x2 i ^ , -56x^-64x-8 
'^S - -3{T3^^^2 l^i - -J + 3{l-x)4 in ^ 3{l-x)4 H , 



120X+8 
3(l-x)2 ' 



^6 = 3(1^^^2(1 -i) + -fz^Wx + fg^lnx + 

7 _ 15x3-|-72x2-35x r ■ /i i 68x3+66x2-38x i^2 ^ 

- 3(r:^p ^^2l^i--J H 3(Yz^)5 in X 

_l_ -150x"^-3534x3+561x2+971x-152 !„ ™ _i_ -3857x3-21960x2+12501x-508 
54(l-x)5 "^-^ 324(l-x)4 ' 

7, _ 64x2+48x7- • /i lA _i_ 136x2+152x 1 2 ^ , -20x3 -404x2 +26x+14 i^ ^ , -127x2-220x-37 
'^S - 3(l-x)3 ^^2 - 3(l-x)4 in 3: i 3(l-a;)4 in 2: i 3(l-x)3 ■ 

Appendix C 

This appendix is devoted to presenting the SSM matching contributions originating from 
the quartic squark vertex proportional to the strong couphng constant ag'- 



92KtsKtb A,B,C=1 7'=1 "^x7 ^ V /^o 

2 

2 

-?3(^A/, Zci) + -g4(^A/, ^c/) 



x{K');K') 



C3 



3 



(106) 
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E E 

A,B,C=1 1=1 



2 M2 



4 



X 



— ^ P° m^- P° 



X 







BC 



In- 



1 X 



X 



8^1 



Ml 



,,2 



X 



QgiK^Ktb mj a,b,c=i 



1 X 



(109) 



In 



1 X 



X 



A2 



r 



C3 



C3 



gi('y^,^^c) - gg2(^^A,^^c) 



(110) 



where = r^^F^^t _ r^^r^^t and P° = T^^F^^t _ r^flp^^t. The mass ratios are 



denoted as before: Zai — m\ /mi-, WAi — rn^ /ttiIo and VA — m^ /m^. 



The exphcit expressions for the functions qi{x, y) read 



(l2(x,y) 
Q3{x,y) 
(l4(x,y) 



2 


x'^ Inx 


3/2 Iny 


+ 


2x'^y'^+5x^y+5xy'^-x'^-y'^-22xy+5x+5y+2 


(111) 


3{x-y) 


[{1-x)^ 


(1-J/)*J 




2 


xlnx 


2/ In J/ 


+ 


-x'^y'^+5x'^y+5xy^+2x'^+2y^-10xy-7x-7y+ll 


(112) 


3{x-y) 


[{1-xr 


(1-J/)^J 


9(l-j;)3(l-2/)3 ' 


4 


x^ Ins 


2/2 In y 


+ 


-6xy+2x+2y+2 


(113) 


3(a;-2/) 


[{1-x)^ 


(l-J/)3j 


3(l-x)2(l-j/)2 ' 


4 


xlnx 




+ 


-2xy-2x-2y+6 


(114) 


3{x-y) 






3(l-x)2(l-2/)2 • 
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